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§ ' Summary 

^ : 

I In the harmonic description of general relativity, the principle part of Einstein equa- 

■ tions reduces to a constrained system of 10 curved space wave equations for the compo- 
I ' nents of the space-time metric. We use the pseudo-differential theory of systems which are 

in ! 

^ , strongly well-posed in the generalized sense to establish the well-posedness of constraint 
CN ■ 

■ preserving boundary conditions for this system when treated in second order differential 

^ I form. The boundary conditions are of a generalized Sommerfeld type that is benevolent 

O . for numerical calculation. 
CT; 

U '. 

W). 1. Introduction 



Consider a first order symmetric system of partial differential equations 

ut = Aux + Buy -f F 

on the half-space x > 0, — oo < y < oo, t > 0. Here w is a vector valued function with n 
components and A = A*, B = B* are symmetric matrices which depend smoothly on x, y 
and t. Also, A is not singular at x = 0. At t = 0, we give the initial condition 

u{x,y,0) = f{x) 

and, at a; = 0, boundary conditions 

Bu = g 

which are strictly maximally dissipative. In this case one can use integration by parts 
to derive an energy estimate. The result can easily be extended to prove local existence 
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of solutions of quasilinear systems because integration by parts allows us to estimate the 
derivatives. 

This result depends heavily on integration by parts. If the boundary conditions are 
not maximally dissipative or the system is not symmetric hyperbolic, new techniques 
are needed. There is a rather comprehensive theory based on the principle of frozen 
coefl&cients, Fourier and Laplace transform and the theory of pseudo-differential operators 
which give necessary and sufficient conditions for well-posedness in the generalized sense 
[1],[2],[3],[4],[5, Chap. 8], [6, Chap 10], [7]. This theory can also be extended to second order 
systems. In the context of pseudo-differential operators one can - in the same way as for 
ordinary differential equations - write a second order system as a first order system [8] . 

To make the results of the above theory more precise and how to apply it we have 
included an Appendix, which should make it easier to read the paper. 

In this paper we shall demonstrate that flexibility of the permissible boundary condi- 
tions can be applied to solve the constraint problem for the harmonic Einstein equations. 
The pseudo-differential theory is applied here to the second order harmonic formulation, 
which was used to establish the first well-posed Cauchy problem for Einstein's equations 
[9]. 

The importance of a well-posed constraint preserving initial-boundary value problem 
(IB VP) to the simulation of Einstein equations has been recognized in numerous recent 
works. The pseudo-differential treatment of the IBVP presented here is most similar to 
treatments of first order formulations by Stewart [10], Reula and Sarbach [11] and Sarbach 
and Tiglio [12]. A well-posed IBVP for the nonlinear harmonic Einstein equations has been 
formulated for a combination of homogeneous Neumann and Dirichlet boundary conditions 
(or boundary data linearized off these homogeneous conditions) [13]. Here we consider 
strictly dissipative, Sommcrfcld-type boundary conditions for these harmonic equations, 
which have proved to be more robust in numerical tests [14,15] . The only general treatment 
of the nonlinear case has been given by Friedrich and Nagy, based upon a quite different 
first order formulation of the Einstein equations [16] 

In the next section we will explain our theory for the wave equation with boundary 
conditions which cannot be treated by integration by parts. We apply the results in the 
third section to the constraint problem of the linearized Einstein equation. 

In the last section we will give a more physical interpretation of our technique and 
explain how it applies to the full Einstein theory. 
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2. Well-posed problems 

Consider the half-plane problem for the wave equation 

vtt = + Vyy + F, < X < OO, — oo < y < oo, t >0, (1) 

with smooth bounded initial data 

v{x, y, 0) = fi{x, y), Vt{x, y, 0) = f2{x, y), (2) 

and boundary condition 

vt = OLVx + ^Vy + g, a; = 0, — oo < y < oo, t > 0. (3) 

Here a, /3 are real constants. We assume always that a > 0. Also, all data are C°°— smooth 
compatible functions with compact support. 

The usual concept of well-posedness is based on the existence of an energy estimate 
which is often derived by integration by parts. Let F = g = and 

/• oc POO 

\\v{;;t)f= / / \v{x,y,t)\''dxdy 

Jo J -oo 

denote the usual L2-norm. Integration by parts gives us , for the energy E, 

d d 

-E =: -Q-^{\\vt\? + + \\vyf) = -2 y vt{Q, y, t)vxiO, y, t)dy 

If in (3) /3 7^ 0, then there is no obvious way to estimate the boundary flux in terms of 
E. Of course, if = and a > 0, there is an energy estimate. To be able to discuss the 
general case we will define well-posedness as in the Appendix. 
We start with a very simple observation: 

Lemma 1. The problem is not well-posed if we can find a solution of the homogeneous 
equation (1) which satisfies the homogeneous boundary condition (3) and which is of the 
form 

v{x,y,t) = e'*+''^ycp{x), |(^|oo < oo. (4) 

Here ^p{x) is a smooth bounded function and ui is real and s — r] a complex constant 
with T] — Re s > 0. 

Proof. If (4) is a solution, then 



is also a solution for any 7 > 0. Thus we can find solutions which grow arbitrarily fast 
exponentially. 

We shall now discuss how to determine whether such solutions exist. We introduce 

(4) into (1),(2) and obtain 

(fxx - {s"^ + = 0, (5) 
s(p{0) = aipxiO) + iPu;ip{0), \ip\oo < 00. (6) 

(5) , (6) is an eigenvalue problem. We can phrase Lemma 1 also as 

Lemma 1'. The problem is not well-posed if (5), (6) has an eigenvalue with Re s > 0. 

(5) is an ordinary differential equation with constant coefficients and its general solution 
is of the form 

i) = (jie'^i^ + (J2e"^", (7) 
where ki = +y/ + oj^, K2 = — V + uj^ are the solutions of the characteristic equation 

- (s^ + oj^) = 0. 

We define the ^ by 

— TT < arg(s^ + cj^) < TT, arg\/ s'^ + lo"^ = ^arg(s^ + o;^). (8) 

Thus 

RcKi > and ReK;2 < for Re s > 0, respectively. (9) 

By assumption ?) is a bounded function and cannot contain exponentially growing compo- 
nents. Therefore, by (9), this is only possible if ai = 0, i.e. 

= (726'^^^ (10) 

Introducing (10) into the boundary condition gives us 

(s — aK2 — ij3uS)a2 = 0. (11) 

Re s > and ReK;2 < tell us that there are no solutions for Re s > since, by 
assumption, a > 0. 



- 5- 



For the purpose of proving well-posedness, from now on we assume that the initial 
data (2) vanish. This may always be achieved by the transformation 

u = v-e-'h-te-\f2 + h). 

(We start the time evolution from 'rest'.) Then we can solve (l)-(3) by Fourier transform 
with respect to y and Laplace transform with respect to t and obtain the inhomogeneous 
versions of (5) and (6), 

Uxx — (s^ + uj^)u = —F, u = u{x, uj, s), 

(12) 

{s — aK,2 — i(3u!)u{0, u), s) — g{u), s). 

Since (5) and (6) have no eigenvalues for Re s > 0, we can solve (12). By inverting the 
Laplace and Fourier transform, we obtain a unique solution. 

It is particularly simple to calculate the solution for F = 0. Corresponding to (10) 
and (11) we obtain 

ii = e'^2^-fi(0,a;,s) (13) 

where 

(s — aK2 — iPu>)u{0, u>, s) = g{u>, s). 
To obtain sharp estimates we need two lemmas. 
Lemma 2. There is a constant 6i > such that 



Proof. 



Let 



Then 



ReK = Re^/oJ^~-\-~s^ > diTj, s — i^-\-r], Re s = rj. (14) 



' = VW-(e)^ + 2i^V + (a;')' + (O' = 1- 



If T]' » 1, then k' Ri rj' and (14) holds. Thus we can assume that {u')^ + {^)^ + {rj')^ < 
const. Assume that (14) is not true. Then there is a sequence 

^' C Co, v' Vo 
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such that 



Re — ^ 0. (15) 



This can only be true if r/^ = 0. If (a;^)^ > then (15) cannot hold. If {u'^f < (Co)^ 

then (^')^ > I and, for sufficiently small r]', (8) gives 



Co 



(16) 



and (14) holds. The same is true if (Co) = (ooq) ■ This proves the lemma. 

Lemma 3. Assume that a > and |/?| < 1. There is a constant ^2 > such that, for all 
uj and s with Re s > 0, 

\s — aK2 — i(3u)\ > S2\^\s\'^ + (17) 
We use the same normalization as in Lemma 2 and write (17) as 



\L\ =: - I3u') + V + a^J{u'Y-{^'Y + 2^i'r]' + {r]'Y\ 
><^vW + W + W, (a;')' + (CO' = 1- 



(18) 



Since a > and Re Vs^ + a;^ > 0, the inequality holds for r]' » 1. Thus we can assume 
that \r]'\ < const. Assume now that there is no 5 > such that (18) holds. Then there is 
a sequence Co' — <^0' ^' ^ such that 



(19) 



Using (8) we obtain 

^{i'o - + c^vw - n'o)' if > eo 



Clearly, Lq 7^ if < 1. Thus we arrive at a contradiction and (17) holds. This proves 
the lemma. 

We can now prove 
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Theorem 1. There is a constant K such that the solution (13) satisfies the estimates 

|mx(0,cu,s)| < K\g{uj,s)\, 
y/\s\^+u;^-\u{0,u,s)\ < K\g{io,s)\. 

Therefore we can use the theory of pseudo-differential operators to obtain the estimates 
and results of the Appendix. In particular, the problem is strongly well-posed in the 
generalized sense. 

Proof. By (13) and (17), 

\ux{0,uj, s)\ < \k2\ \u{0,uj,s)\ — I + s^l \u{0,uj,s)\ < K\g{uj,s)\. 

The estimates for the other derivatives follow directly from (17) and (14). 

(A18) is the first order version of (12). For F = 0, we have Ux = y^\s^~l-oj^ v and 
therefore, for a; = 0, 

•\/|sp + \v{0,uj,s)\ = \ux{0,u>,s)\ < K\g{u!,s)\. 

Thus the required estimate (A20) holds and we can apply the Main theorems Al and A2 
and obtain the estimates (A21) and (A22). This proves the theorem. 

3. Linearized Einstein equations 

We consider the half-plane problem for the linearized Einstein equations 



i-d^ + dl + dl + dl) 















j^y 




ryty 




r^yy 


r^yz 











(26) 



a; > 0, t > 0, — oo < y < oo, — oo < z < oo, 
together with the constraints C", 

c* = dtY' + dxi'"' + dy^'-y + a,7'' = 0, 
= 5^7*" + a,7"" + dy^'^y + a,7"^ = 0, 
cy = da*y + 9^7^^ + dy^yy + d,^y' = 0, 

= da^' + a^7^^ + dyY^ + d^Y^ = 0. 

The constraints are also solutions of the wave equation. We can guarantee that they remain 
zero at later times if = 0, a = {t, x, y, z), are part of the boundary conditions for (26) 
at a; = 0. 



(27) 
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A possibility is 
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(28) 



Here ai, 02, as, bi, 62 and ci, C2 are real constants such that the eigenvalues of 




1 
bi 62 



and 



1 

Cl C2 



(29) 



are real and negative. 
We want to prove 

Theorem 3. The half-plane problem for the system (26) with boundary conditions (28) 
is well-posed in the generalized sense if the eigenvalues of the matrix (29) are real and 
negative. 
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Proof. We can assume that F = and need only show that the estimate of Theorem 1 
holds for every component 7'^'^ with |^| denoting the Euclidean norm of all components of 
the forcing. 

Fourier transform the problem with respect to the tangential variables and Laplace 
transform it with respect to time. For every component 7*-' we obtain 

{dl - (s^ + + ujD) = 0, I7I00 < 00, (30) 

which are coupled through the corresponding transformed boundary condition. We start 
with the last component 

(a^-(.2 + a;? + a;i))r^ = 0, a; > 0, 

(ol) 

= (^^^^^ + for x = 0, |7^^|oo < 00. 

This is a problem which we have treated in the last section. By Theorem 1, we gain one 
derivative on the boundary. For 7^'^, 7^^ we have the same result. 
The coupled boundary conditions for 7*^, 7^^ 

(0 l\a Ci^^^i.. (r'^^f a" 



can be decoupled. There is a unitary matrix U such that 



Introducing new variables by 



gives us 

^ if) - (0 th CoO CO - ip 

We start with the last component. By (31), we have gained a derivative which we loose by 
calculating 9^7^^ = iu>2^^^- But then we gain a derivative by solving for 7^^. The same is 
true for 7*^. The process can be contimied and the theorem follows. 

Clearly, Theorem 3 is also valid when the matrices of (28) for the tangential derivatives 
are strictly upper triangular, i.e., only terms above the diagonal are not zero. This allows 
the sequential argument associated with (31) and (32) above. One can also generalize the 
result to full matrices provided the elements are sufficiently small. 
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4. Constraint preserving boundary conditions in Sommerfeld form 

The example in the preceding section illustrates how the pseudo-differential theory 
can used to establish a constraint-preserving IBVP for the linearized Einstein equations 
which is well-posed in the generalized sense. There are further boundary conditions not 
covered in the example that can be treated by the same technique because the boundary 
conditions can be written in the form (28) with the matrices for the tangential derivatives 
strictly upper triangular. Here we consider a simple hierarchy of Sommerfeld boundary 
conditions. The geometric nature of the construction is more transparent using standard 
tensor notation based upon spacetime coordinates a;" = {t,x,y,z). Our results center 
about systems whose components satisfy the scalar wave equation, now written in terms 
of the Minkowski metric 77*^^ = diag{—l, 1, 1, 1) as 

V^'^dadf,^ = {-dl + dl + dl + dl)^ = 0, 

in the half-space x > 0, t > 0, —00 < y < 00, —00 < z < 00. 
For this IBVP, the energy 

E=IJ (^{dt^f + {d,^f + {dy^f + {d.^f^xdydz 

satisfies 

dtE = - / J^dydz 

Jx=0 

where the energy flux through the boundary is 

This leads to a range of homogeneous, dissipative boundary conditions 

Adt + Bda^^ = 0, 

subject to values of A and B such that > 0. The Dirichlet boundary condition 
[A = 1,B = 0) corresponds to the case where $ has an odd parity local reflection symmetry 
across the boundary; and the Neumann boundary condition {A = 0, B = 1) corresponds to 
an even parity local reflection symmetry. Both the homogeneous Dirichlet and Neumann 
conditions are borderline dissipative cases for which = 0. The description of a traveling 
wave carrying energy across the boundary requires an inhomogeneous form of the Dirichlet 
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or Neumann boundary condition, which provides the proper boundary data for the wave 
to pass through the boundary. However, in numerical simulations, such inhomogeneous 
Dirichlet or Neumann boundary data can only be prescribed for the signal and the numer- 
ical error is reflected by the boundary and accumulates in the grid. This can lead to poor 
performance in simulations of dynamical nonlinear systems such as Einstein's equations. 
For such computational purposes, it is more advantageous to use the strictly dissipative 
Sommerfeld condition given hy {A — 1, B — —1), i.e. 

{dt - d,)^ = 0, 

for which = (5^$)^. The Sommerfeld condition is based upon the characteristic direction 
in the 2-space picked out by the outward normal to the boundary and a timelike direction 
(the evolution direction) tangent to the boundary. 

Before formulating Sommerfeld boundary conditions for the constrained harmonic 
Einstein equations, it is instructive to consider the analogous case of Maxwell's equations 
for the electromagnetic field expressed in terms of a vector potential A^^ = (A*, A^^A^^A^) 
(see e.g. [17]). Subject to the Lorentz gauge condition 

C := d^A^' = 0, (33) 

the Maxwell equations reduce to the wave equations 

v'^^dadpA'^ = 0. (34) 

The constraint C also satisfies the wave equation ij^f^dadpC = 0. Thus Cauchy data 
A^\t=o and dtA^\t=o which satisfies C\t=o = dtC\t=o = 0, i.e. for which C also has vanishing 
Cauchy data, leads to the well-known result that the constraint is preserved for the Cauchy 
problem. 

In order to extend constraint preservation to the IBVP, the boundary condition for 
A'^ must imply a homogeneous boundary condition for C. One way to accomplish this is to 
use locally reflection symmetric boundary data, as in the scalar case. For instance, for the 
above half-space problem, the even parity boundary conditions dxA^\x=o = dxAy\x=o = 
dxA^\x=o = A^lx^o = imply that dxC\x^o = 0. These boundary conditions are dissipa- 
tive and lead to a well-posed IBVP for Maxwell's equations. The analogous approach has 
been used to formulate a well-posed constraint-preserving IBVP for the nonlinear Einstein 
equations [13]. However, this leads to a combination of Dirichlet and Neumann boundary 
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conditions which is only borderhne dissipative and the results of tests for the nonlinear 
Einstein problem show that a strictly dissipative Sommerfeld-type boundary condition 
gives better numerical accuracy [15]. 

The pseudo-differential theory offers alternative approaches. Consider the IB VP with 
Sommerfeld boundary conditions 

(dt-d,)iA' + A-)=g, (35) 

(dt - d,)Ay = gy, (36) 

{dt - d,)A^ = g\ (37) 

^{dt - d,){A' - A^) + dt{A* + A-) + dyAy + d,A' = ^g, (38) 

where g, gy and g^ are free Sommerfeld data. The IBVP is well-posed in the generalized 
sense if the estimates of Theorem 1 are satisfied. Using the argument associated with 
(31) in the preceding section, it follows from (35)-(37) that A* + A^, A^, A^ and their 
derivatives satisfy these estimates. Next, consider the Sommerfeld boundary condition 
(38) for A'^ — A^. Using the argument associated with (32), the tangential derivatives 
of {A^ + A^), Ay and A^ introduce no problem, and the estimates extend to A* — A^ 
and its derivative. Thus the Sommerfeld boundary conditions (35)-(38) guarantee a well- 
posed IBVP for the system (34). But (35)-(38) also imply that the constraint satisfies the 
homogeneous boundary condition C = 0, as is evident by rewriting (33) in the form 

C=lidt- d,)iA' - A^) + l{dt + d,)iA' + A^) + dyAy + d,A\ (39) 

Although the boundary conditions (35)-(38) lead to a well-posed, constraint pre- 
serving IBVP for Maxwell's equations, they do not correspond to any physically familiar 
boundary conditions on the electric and magnetic field vectors E and B. However, there 
are numerous options in constructing boundary conditions by this approach. For instance, 
consider the choice 

{dt-d^){A' + A-) = Q, (40) 
{dt-d,)Ay + dy{A' + A-) = Q, (41) 
{dt - + dz{A' + A^) = 0. (42) 

]^{dt - d^){A' - A^) + dt{A' + A^) + dyAy + d,A^ = 0. (43) 



Again, the Sommerfeld condition (40) implies that A* + = and its derivatives satisfy 
the estimates required for Theorem 1. In addition, the sequential manner in which the 
tangential derivatives of previously estimates quantities enter (41)-(43) ensures that all 
components of A^^ satisfy the required estimates. Furthermore, (40)-(43) imply that the 
constraint satisfies the homogeneous boundary condition C = 0, as is again evident from 
(39). By using the wave equations (34) and the constraint (33), it can be verified that 
these boundary conditions give rise to the familiar plane wave relations = —B^ and 
— on the components of the electric and magnetic fields. These relations imply 
that the Poynting flux E x B leads to a loss of electromagnetic energy from the system. 

The half-space problem for the linearized harmonic Einstein equations can be treated 
in a similar manner. Given a metric tensor g^^ with inverse g^^ and determinant g, the 
linearized Einstein equations imply that the perturbations = S{^/—gg'^^) satisfy the 
10 wave equations 

ri'^^dadp-i'''' = (44) 

provided the 4 constraints 

C^" a^7^^ = (45) 

are satisfied [17]. The constraints constitute the harmonic gauge conditions which reduce 
the Einstein equations to a symmetric hyperbolic system of wave equations. 

A simple formulation of Sommerfeld boundary conditions for the half-space problem 
< a; < oo for the linearized gravitational field can be patterned after the above treatment 
of the electromagnetic problem. For convenience we write = {y,z). First we require 
the 6 Sommerfeld boundary conditions 

{dt - d,) = q^^, (46) 

(5t-a,)(y^ + 7"^) = 9^ (47) 

(at-c».)(7" + 27'" + 7"")=9, (48) 

where q"^^ , q^ and q are free Sommerfeld data. Then the constraints are used to supply 
4 additional Sommerfeld-type boundary conditions in the hierarchical order 

= lidt- a,)(7*^ - + dt{-f'^ + + dBl^'' - ig^^ = 0, (49) 
+ = ^{dt - a,)(7« - 7^-) + at(7« + 27*- + 7^-) + 83(1''' + t""") -^9 = 0' (50) 
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= \{dt- + in + 9^(7" + f + dBl''' - = 0. (51) 

In expressing the constraints in this form, we have used (47)-(49) and the prior constraints 
in the hierarchy. 

The sequence of Sommerfeld conditions (46)-(51) for [dt — dx)l'^^ have the property 
that all tangential derivatives of 7^*^ only involve prior components in the sequence. This 
allows us to again use the arguments associated with (31) and (32) to obtain the estimates 
for all components of 7'^'^ and their derivatives which are required for Theorem 1. 

The boundary conditions (46)-(51) offer a simple and attractive scheme for numer- 
ical use. The pseudo-differential theory allows many more possibilities for a well-posed, 
constraint-preserving IBVP. It would be of value to find a version with a simple phys- 
ical interpretation in the analogous way that (40)-(43) implies a positive Poynting flux 
in the electromagnetic case. However, this issue is complicated by the lack of a unique 
expression for the gravitational energy flux except in the asymptotic limit of null inflnity. 
One practical alternative for a numerical scheme would be to use an external solution to 
provide the Sommerfeld data on the right hand sides of (46)-(48). This data can be ob- 
tained either by matching to an external linearized solution or, in the nonlinear case, by 
Cauchy-characteristic matching [18]. 

Our results generalize to the curved space linearize harmonic wave equation, whose 
principle part has the form 

g'^^do^dp^'', (52) 

determined by a given space-time metric g'^^{t^ x, z), i.e. a matrix which can be trans- 
formed at any point to diagonal Minkowski form. This variable coefficient problem is 

well-posed in the generalized sense if all frozen coefficient problems are well-posed [5]. 
This result is insensitive to lower order terms in the equations. This principle of frozen co- 
efficients is an important result of the pseudo-differential theory. We can reduce the frozen 
coefficient problem to the above Minkowski space problem by adapting the harmonic co- 
ordinates so that the boundary is given by a; = and then introducing a orthonormal 
tetrad, 

where X^ is in the direction Va^c normal to the boundary. Freezing the tetrad at a 
boundary point, we can then introduce the linear coordinate transformation 

t — TqiX , X — XqiX , y — Yqi^x , z — Zq^x . 
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In the 5" coordinates, the frozen coefficient problem reduces to the Minkowski space 
problem, which we have treated. 

The full treatment of the Einstein equations requires taking into account the relation 
1^^ = V~99^'^-> which converts (52) into a quasi-linear operator. In this case the pseudo- 
differential theory outlined in the Appendix establishes that the IBVP is well-posed locally 
in time. 

Appendix 

In this section we shall give a short summary of the theory for first order systems. 
This includes also second order systems because our theory is based on pseudo-differential 
operators and therefore one can always write a second order system of differential equations 
in terms of a first order system of pseudo-differential operators. (See the references in the 
Introduction.) 

Consider a first order system 

m 

ut = P{d/dx)u + F, P{d/dx) = A d/dxi + ^ Bjd/dxj {Al) 

with constant coefficients on the half-space 

t >0, xi >0, — oo < Xj < oo, j = 2, . . . ,m. 

Here u{x,t) = (u^^\x,t), . . . ,u^'^\x,t)) is a vector valued function of the real variables 
(x, t) = (xi, . . . , Xm, t) and A, Bj are constant n x n matrices. 

We assume that the system is strictly hyperbolic, i.e., for all real oj = (a;i,a;_), a;_ = 
{u!2, ■ ■ ■ ,(^m) with I a; I = 1, the eigenvalues of the symbol 

m 

P{iuj) = iAuJi +iB{uj-), S(cj_) = ^ Bjujj, {A2) 

are purely imaginary and distinct. 

In [1] the theory has been extended to the case where the eigenvalues have constant 
multiplicity and there is a complete set of eigenvectors. In particular, the theory holds if 
the system consists of strictly hyperbolic subsystems which are only coupled through lower 
order terms and the boundary conditions. 
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We assume also that A is nonsingular and without restriction we can assume that it 
has the form 

--{-^' ;-)• 

Here A^, A^^ are real positive definite diagonal matrices of order r and n — r, respectively. 
For the singular case, sec [7]. 

For t = 0, we give initial data 

uix,0) = f{x) (^4) 

and for, xi — 0, r boundary conditions 

u\0,x-,t) = Su^\0,x-,t) + g{x-,t), X- = {x2, . . . .iXm)- (Ad) 

All data are smooth, compatible and have compact support. 

The usual theory for well-posed problems depends strongly on the assumption that 
the system is symmetric and that the boundary conditions are maximally dissipativc. If 
any of those two arguments is not satisfied, then the theory docs not give anything. We 
want to discuss a concept of well-posedness for which we obtain necessary and sufficient 
conditions. 

The main ingredient of a definition for well-posedness is the estimate of the solution 
in terms of the data. (See [5, Sec. 7. 3]). We will consider (Al),(A4),(A5) with homogeneous 
initial data / = and use 

Definition Al. Let f{x) = 0. We call the problem strongly well-posed in the generalized 
sense if, for all smooth compatible data F, g, there is a unique solution u and in every time 
interval <t <T there is a constant Kt which does not depend on F and g such that 

f\\ui;T)rdT+ f\\ui;T)\\ldT<KT{f\\Fi.,T)fdT+ ^ t) |p_dr}. (^6) 

Jo Jo Jo Jo 

Here || ■ ||, || • ||_ denote the L2-norm with respect to the half-space and the boundary space, 

respectively. 

We start with a simple observation. For F = = 0, we construct simple wave solutions 

m 

u{xi,x-,t) ^e'*+'^'^'''>-ifi{xi), (cj,a;)_ = ^cjjXj-, {A7) 

satisfying the boundary conditions 

^■^(0, X-) = S(p^\0, X-), \(p\oo < oo. (AS) 

We have 



Lemma Al. The half-plane problem is not well-posed if, for some ojq and complex sq 
with Re So > 0, there is a solution (A7) which satisfies (A8). 

Proof. If there is a solution then, by homogeneity, 

is also a solution for any 7 > 0. Thus there are solutions which grow arbitrarily fast 
exponentially. This proves the lemma. 

We shall now derive algebraic conditions so that we can decide whether these simple 
wave solutions exist. 

Introducing (A7) into (Al) and (A5), we obtain 

Sip = Aifx + iB(uj-)(p, a; > 0, 

{A9) 

ip\0) = V^(0), l^loo < 00. 

(A9) is an eigenvalue problem for a system of ordinary differential equations which can be 
solved in the usual way. Let k, denote the solutions of the characteristic equation 

Bet\AK- (si -iB(u;-))\ = 0. (AlO) 

One can prove 
Lemma A2. 

1 ) For Re s > 0, there are no k with Re k = 0. 

2) There are exactly r eigenvalues with Re k < and n — r eigenvalues with Re k > 0. 

3) There is a constant 5 > such that, for all s = i^ + r], ^, 77 > 0, and all 

|Re k\ > Sr], s — + rj, ^, ry > real. 

(See [2]). 

We can now write down the general solution of (A9). If all kj are distinct, the solution 
is of the form 

E (^je"'"hj+ Yl <'3^"'"hj. (All) 

Here hj are the corresponding eigenvectors. (If the eigenvalues are not distinct, the usual 
modifications apply.) 
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Since we are only interested in bounded solutions, all aj in the second term are zero. 
Introducing (p into the boundary conditions at a; = gives us a linear system of r equations 
for r unknowns (cri, . . . , cr^) = a which we write as 

C(oj_,s)a = 0. (A12) 

The problem is not well-posed if for some a;_, there is an eigenvalue sq with Re sq > 0, 
i.e., DetC(a;_, sq) = 0. Then the linear system of equations (A12) and therefore also (A9) 
has a nontrivial solution. 

From now on we shall assume that Det C 7^ for Re s > 0. Then we can solve the 
initial boundary value problem by Laplace transform in time and Fourier transform in 
the tangential variables. For convenience, we start the solution from 'rest', i.e., u{x,0) = 
f{x) = 0. Then we obtain 

su = Aux + iB(u>-)u + F, 

(A13) 

Since, by assumption, (A9) has only the trivial solution for Re s > and |Re«:| > drj, 
(A12) has a unique solution. Inverting the Fourier and Laplace transforms we obtain the 
solution in physical space. 

It is particularly simple to solve (A13) if F = 0. By (All) and (A12), 

■u(0,a;_,s)= ^ Ujc'^'^hj, 

Re Kj <0 

where the cr^ are determined by 

C{uj-,s)a = g. 

Definition A2. Consider (A13) with F = 0. The problem is caUcd boundary stable if, 
for all Lv,s with rj = Re s > 0, there is a constant K which does not depend on u>,s and g 
such that 

\u{0,u},s)\ < K\g{oj,s)\. (AU) 

One can also phrase the condition as: The eigenvalue problem (A9) has no eigenvalues 
for Re s > or Det C(a;_, s) for Re s > 0. 

The estimate (A14) is crucial to the theory. It allows us to construct a symmetrizer 
to obtain an energy estimate in the generalized sense for the full problem. 
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We introduce normalized variables 

s' = s/Vl^P + li^-P = 'i-C + v'l ui'_ = U-/^/\s^~+]uJI\^ 
and write (A13) as 

-Aua; + a/|s|2 + |a;_|2 (s'l - iB(u;'_)) u = F, 

{A15) 

u^{0) - Su^^iO) = g. 

Main theorem Al. Assume that the half-plane problem is boundary stable. Then there 
exists a symmetrizer R = R{s',(jj'_) with the following properties. 

1) R is a smooth bounded function of s', uj'_ and the coefficients of A, Bj and S. 

2) RA is Hermitian and for all vectors y which satisfy the boundary conditions 

{y,RAy)>6i\y\''-c\g\\ 

3) y|s|2 + |a;_|2Re{A {s'l - iB{uj'_))} > 52r]I. 

Here 5i, ^2 > 0, c > are constants independent of s', uj'_. 
We can now prove 

Main theorem A2. Assume that the half-plane problem is boundary stable. Then it is 
well-posed in the sense of Definition Al. 

Proof. Multiplying (A15) by R we obtain, for the scalar product with respect to xi, 

Re{u, RF) ^Re{-{u,RAdu/dxi) + (u, + |a;_|2 • R{s'I - iB{iu'_)uj} 

= Bje{~{u, RAu)\2^^ + (u, RisI - zS(a;_)«) } 
> 511^(0, s,U-)\'^ + S2r]\\u{xi, - c\gf. 

Thus we obtain 

r]\\u{xi,s,u;)f + \u{0,s,u;)\'^ < const. + c|^p^ . (AlQ) 

Inverting the Laplace and Fourier transform proves the theorem. 

Remark. The estimate (A14) and the properties of the symmetrizer and therefore also 
the estimate (A16) need only be valid for ?7 > ?7o > 0. This is important if lower order 
terms or variable coefficients are present. 
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One could have derived the estimate (A16) by directly calculating the solution of 



and (A9) as an equation for pseudo-differential operators. The theory of pseudo-differential 
operators has far reaching consequences. In particular, the computational rules for pseudo- 
differential operators show: 

1) The estimate is also valid if the symbols depend smoothly on x,t, provided we assume 
that 1] > r]o, rjo sufficiently large. Here rjo depends on a finite number of x, t-derivatives of 
the symbols. Therefore we extend the estimate to systems with variable coefficients. 

2) Well-posedness will not be destroyed by lower order terms. Therefore one can localize 
the problem, and well-posedness in general domains can be reduced to the study of the 
Cauchy problem and the half-plane problems. 

3) The principle of frozen coefficients holds. 

4) The properties of the pseudo-difi^erential operators allow us to estimate derivatives in 
the same way as for standard partial difi^erential equations. Therefore we obtain well- 
posedness in the generalized sense for linear systems with variable coefficients which gives 
us the corresponding local results for quasi-linear problems. 

Since pseudo-differential operators arc much more fiexible than standard differential 
operators, we can always write second order systems as first order systems. Consider, for 
example, the problem we discuss in Section 2. 



(A9). However, the importance of the symmetrizer is that we can consider ^ as a symbol 



utt = Uxx + Uyy + F for X > 0, — oo < y < oo, t > 0, 
Uf — aux — Puy = g, X = 0, — oo < y < oo, t >0. 



{A17) 



After Laplace-Fourier transform it becomes 



sit — aux — Piuju = g. 



Introducing a new variable itx = 



•\/|sP -\- uj'^v gives us the first order system 




(A18) 



s'u — av — j3iuj'u = g, 
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where 

' 



Definition A3. We caii tie problem (A17) strongly well-posed in the generalized sense 
if the corresponding Grst order problem (A18) with general data F, g has this property. 

If (A18) is boundary stable, i.e., if the estimate (A14) holds with g replaced by g, then 
we can use the same technique as in [2, Sec. 4] to construct the symmetrizer. Therefore 
Main Theorems 1 and 2 are valid and we obtain the estimate (A16) with F, g replaced by 

Starting from (A17), F,g satisfy (A19) and therefore (A14) becomes 

\u{0,s,u)\ + \v{0,s,uj)\ < f |^(a;,g)|. (^20) 

In Section 2 we prove that (A20) holds by directly calculating the solution of (A17) for 
F = 0. 

Similarly, in terms of F and g for the second order system, the estimate (A16) becomes 
{\s\^ + u^){r]\\u{-,LO,s)f + \u{0,s,Lo)\^)< const. (^^||Ff + cj^p^ . {A21) 
We have also 

\\u^f = ilsf + u^) \\vf < {\sf + u;^)\\uf. {A22) 

Therefore, by inverting the Laplace and Fourier transform, we can estimate the L2-norm 
of all first derivatives in terms of the L2-norm of the data. Thus we gain one derivative. 
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